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Classical chaotic systems are distinguished by their sensitive dependence on initial 
conditions. The absence of this property in quantum systems has lead to a number 
of proposals for perturbation-based characterizations of quantum chaos, including 
linear growth of entropy, exponential decay of fidelity, and hypersensitivity to per- 
turbation. All of these accurately predict chaos in the classical limit, but it is not 
clear that they behave the same far from the classical realm. We investigate the 
dynamics of a family of quantizations of the baker's map, which range from a highly 
entangling unitary transformation to an essentially trivial shift map. Linear entropy 
growth and fidelity decay are exhibited by this entire family of maps, but hypersen- 
sitivity distinguishes between the simple dynamics of the trivial shift map and the 
more complicated dynamics of the other quantizations. This conclusion is supported 
by an analytical argument for short times and numerical evidence at later times. 



I. INTRODUCTION 

A full characterization of quantum chaos is an elusive matter. Classical chaotic sys- 
tems are distinguished by their exponential sensitivity to initial conditions. Quantified in 
terms of Lyapunov exponents, this characterization is the key ingredient in any definition 
of classical chaos. The linearity of quantum mechanics prohibits such sensitivity to initial 
conditions, thus obstructing any straightforward extension of the classical definition of chaos 
to quantum systems. The standard fix is to categorize as "chaotic" those quantum systems 
that are chaotic in a classical limit. These systems are not strictly chaotic by the classical 
definition — they are quasiperiodic — but they have properties, involving the spectrum of en- 
ergy eigenvalues and the behavior of energy eigenstates, that are distinctly different from 
that for quantizations of classically regular systems. 
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So far, however, little agreement has been reached on a characterization of quantum chaos 
that does not make reference to a classical limit. The classical approach of looking at small 
perturbations of the initial state fails due to the unitarity of quantum dynamics. Attempts 
at a direct dynamical characterization of quantum chaos, which applies even in the hard 
quantum regime, thus generally look at the effects of small perturbations of the dynamics. A 
number of perturbation-based criteria (or signatures) have beenproposed, including linear 
growth of entropy in the presence of stochastic perturbations |l[ |2|; exponential decay of 
fidelity between states that evolve under two close, but distinct unitary transformations 0,1^; 
and hypersensitivity to perturbation, which considers stochastic perturbations and compares 
the amount of information known about the perturbation to the resulting reduction in the 
system entropy 0, IE 013; a relation called the information- entropy trade-off. 

Typically, a quantum system whose classical limit is chaotic exhibits all three of these 
signatures. They are clearly inequivalent, however. For instance, hypersensitivity to pertur- 
bation can be viewed as a measure of how fast, how widely, and how randomly the set of 
all possible perturbed states spreads through Hilbert space. The other two criteria, though 
they report on how widely the perturbed vectors are dispersed through Hilbert space, are 
not sensitive to the exact way in which Hilbert space is explored by the perturbed dynam- 
ics. It is therefore conceivable that there are quantum systems that are chaotic with respect 
to one criterion, but regular with respect to another. In this paper we compare the three 
perturbation-based criteria for a family of quantizations of the quantum baker's map. 

The classical baker's map 0] is a well-known toy mapping whose study has led to many 
insights in the field of classical chaos by demonstrating essential features of nonlinear dy- 
namics. It maps the unit square, which can be thought of as a toroidal phase space, onto 
itself in an area-preserving way. Interest in the baker's map stems from its straightforward 
formulation in terms of a Bernoulli shift on binary sequences. It seems natural to consider a 
quantum version of the baker's map for the investigation of quantum chaos. There is, how- 
ever, no unique procedure for quantizing a classical map; hence, different quantum maps 
correspond to the same classical baker's transformation in the classical limit. The family of 
quantizations used in the present paper is based on the 2 iV -dimensional Hilbert space 
of N qubits. This qubit structure provides a connection to the binary representation of the 
classical baker's map. 

The paper is organized as follows. In Sec. Owe introduce our mathematical notation and 
give precise definitions of the three perturbation-based chaos criteria. Section 11111 reviews 
a family of quantizations for the quantum baker's map. These range from a highly entan- 
gling unitary transformation to an essentially trivial shift map. In Sec. IIVI we give simple 
analytical results concerning the three criteria for the trivial shift map; these results show 
that the trivial shift map exhibits linear entropy growth and exponential fidelity decay and 
suggest that it does not display hypersensitivity to perturbation. Section IVl presents numer- 
ical calculations for the entire family of baker's maps. These calculations show that, unlike 
the other two criteria, hypersensitivity to perturbation differentiates between the different 
quantizations. In Appendix El we formulate a simple model for the form of the information- 
entropy trade-off in the case of vectors distributed randomly on Hilbert space. The model 
serves as a foil for interpreting the results of our numerical work on hypersensitivity. Ap- 
pendix [B] derives the von Neumann entropy of an ensemble of vectors that populate half 
of Hilbert space uniformly; this result is used to bound the information-entropy trade-off 
in the case that the amount of information about the perturbation is one bit. Finally, in 
Sec. EUwe discuss our results. 
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II. CRITERIA FOR QUANTUM CHAOS 
A. Hypersensitivity to perturbation 

1. Definition of hypersensitivity 

In the most general setting, hypersensitivity to perturbation can be defined as follows [11]. 
Consider a system with Hilbert space S, evolving under some unitary evolution and, in 
addition, interacting with an environment with Hilbert space £. Let D and denote the 
Hilbert-space dimensions of the system and environment, respectively. Initially, the joint 
state of the system and environment is assumed to be a product state; i.e., initially there 
is no correlation. After a time t, the joint state of the system and environment is a density 
operator on S®£, which we denote by p to tai- The state of the system at time t, p, is obtained 
by tracing out the environment, 

P = tr £ (p tota i) . (2.1) 
The von Neumann entropy of the system at time t is 

H s = -tr(p\ogp) . (2.2) 

We measure entropy in bits (i.e. we take log = log 2 ). 

Now assume that an arbitrary measurement is performed on the environment. The 
most general measurement 0] is described by a positive-operator- valued measure (POVM), 
{E r }, where the E r are positive operators acting on the environment and satisfying the 
completeness condition 

y~] E r = lg = (environment identity operator). (2-3) 

r 

The probability of obtaining the measurement outcome r is given by 

Pr = tr(/3 total (is® E r )) , (2.4) 

where S.5 is the identity operator on the system. The system state after a measurement 
that yields the outcome r is 

tr e (p total {is ® E r j) 

Pr = —± L ■ (2.5) 

Pr 

We define the system entropy conditional on the outcome r, 

H r = -tr(p r logp r ) , (2.6) 

the average conditional entropy, 

H = Y,PrH r , (2.7) 

r 

and the average entropy decrease due to the measurement, AH = H$ — H. Furthermore, 
we define the average information, 



I = -^PrlogPr 

r 



(2.8) 
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The quantity I is within 1 bit of the minimum average algorithmic information needed to 
specify the measurement outcome r 

Now assume we want to perform a measurement that reduces the average conditional 
system entropy below some given target value, H. We define the quantity 

7 min (#) = inf/, (2.9) 

where the infimum is taken over all POVMs {E r } such that H < H. The function 7 min (if) 
expresses what we call the information- entropy trade-off; it can be interpreted as the min- 
imum information about the perturbing environment needed to keep the average system 
entropy below the target value H . We say the system is hypersensitive to perturbation if 
this information is large compared to the purchased entropy reduction AH = H s — H, i.e., 

» 1 , (2-10) 

in the region of small enough entropy reductions that this ratio reports on the system 
dynamics rather than on the multiplicity of possible perturbations. We characterize this 
region more precisely in the next subsection. 

For the analysis of the present paper, we specialize to the case where initially the system 
is in a pure state, \ipo), an d the unitary system evolution is given by a quantum map B. In 
the absence of any interaction with the environment, the system state after t iterations (or 
time steps) is -B*|^ )- in addition, we assume that the effect of the environment is equivalent 
to a stochastic perturbation. At each time step, a perturbation is chosen randomly from a 
set of unitary maps, {Uk '■ k G /C}, where K, is some index set. 

The joint system-environment density operator after t iterations is then given by 

Ptotal = Pk|^k)(^k| ® P£ , (2.11) 

ke/c* 

where 

|^ k ) = UkMh^B ■ ■ ■ U kl B\ij ) (2.12) 

is the endpoint of a stochastic trajectory labeled by k = (k\, . . . , k t ), pk is the probability 
of k, and the operators P£ are one- dimensional, orthogonal environment projectors. The 
perturbation histories k = (k\, . . . , k t ) are thus recorded in the environment in the form of 
the orthogonal projectors P£. The reduced density operator for the system is given by 

p = tr £ (p tata i) = Pkl^kX^kl • (2.13) 
ke/c* 

In the numerical analysis in Sec. we let the environment be a series of qubits, which 
interact sequentially with the system; hence the perturbation at each time step is a binary 
perturbation, ]C — {0,1}. 

It is in general very difficult to determine the function J min (if) because it is generally 
impossible to do the required optimization over all POVMs. In the numerical results reported 
in Sec. |VJ we restrict the optimization to POVMs of the form 

K= E H, (2-14) 
kei^ 
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where the subsets K r C K} are nonoverlapping subsets of the perturbation histories. Such 
POVMs can be regarded as sampling a coarse-grained version of the perturbation histories. 
For the special form of ptotai considered here, it seems reasonable that ensembles which 
are optimal with respect to this class of measurements are also, to a good approximation, 
optimal with respect to the class of all possible environment measurements. We have not, 
however, been able to prove this statement rigorously. 

The measurements (j2.14j) correspond to forming groups of system vectors iV'k)- Assuming 
that all perturbations are equally likely, i.e., pk = 1/Af for all k, where Af is the number 
of perturbation histories or system vectors the probability of obtaining outcome r is 
given by 

p r = \K r \/Af, (2.15) 

and the system state after a measurement that yields outcome r is the average of the grouped 
vectors, 



tr £ (ptotai(is®£< 



1 



Pr = — * L = hft E iV'kXV'kl • (2.16) 



Pr \K r 



keK r 



We use this simplified framework in the discussion in Sees. IV and V below. 

Even this simplified framework is not sufficient to make the problem tractable for numer- 
ical purposes, because the number of vectors, Af, increases exponentially with the number 
of time steps, rapidly making it impossible to search over all possible ways of grouping the 
vectors l^k)- To get around this, we employ efficient algorithms for grouping the vectors, 
which are plausibly able to find optimal or near-optimal groupings. For the numerical results 
reported in Sec. El we first use a particularly simple, but intuitive grouping algorithm de- 
vised to clarify the procedure and then take a general approach based on genetic algorithms. 
We have compared the results obtained using these algorithms with those obtained using 
other grouping algorithms, some of which have been used previously |7[ , and found that the 
current algorithms are generally superior for the vectors generated by the perturbed baker's 
map. 



2. Quantitative measure of hypersensitivity 

Hypersensitivity to perturbation tests how fast and how fully the state of the perturbed 
system explores the system Hilbert space. To see how this is quantified by the information- 
entropy trade-off, we consider the trade-off relation for vectors that are distributed randomly 
in Hilbert space. Such a relation was formulated in 0, EH, using a model that groups 
the random vectors into spheres of uniform radius (measured by Hilbert-space angle) on 
projective Hilbert space [l4|. We refine this model and its trade-off relation in Appendix [XJ 
The main result is that for Af vectors distributed randomly in d Hilbert-space dimensions, 
the information-entropy trade-off, written in inverse form, is approximated by 

r log Af - I min , log Af > J min > log Af - log d, 

~ | \ogd- -((1 + J min ln2)log(l + J min ln2) - J min ) , log A/" - \ogd > I min > 1. 

(2.17) 

This expression assumes that d is large and that the number of random vectors, though large 
in the sense that Af ^> d, satisfies Af <C 2 d , a situation we refer to as a sparse collection 
of vectors. Examples of the information-entropy trade-off are shown in Fig. |3] (Sec. IV|) . 
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in which the upper solid curve closely resembles the exact situation for random vectors in 
32 dimensions; notice that I m j n follows H in a linear fashion before dropping quickly at 
a "knee" close to the maximum entropy. This agrees with our approximation (|2.17|) for 
random vectors, which is shown as the upper dotted line and the rightmost dotted curve in 
Fig. El 

Several features of the sphere-grouping trade-off (j2.17j) deserve discussion. The number 
of spheres, 2 7min , gives the number of vectors per group Mv — A/"2 _/min . The knee at 7 m i n = 
log M — logd thus corresponds to Mv = d. For J min > log M — logrf, i.e., Mv < d, the 
number of vectors in each sphere is insufficient to explore all the Hilbert-space dimensions. 
This gives a linear dependence on J min , with slope —1 and intercept log A/". In the context 
of a stochastically perturbed map, where M is the number of perturbation histories, this 
part of the trade-off relation tells us about the multiplicity of the perturbation instead 
of about the dynamics of the map. In contrast, for J min < log A/" — logd, i.e., Mv > d, 
where the number of vectors in each sphere is large enough to explore all Hilbert-space 
dimensions, the information-entropy trade-off becomes independent of M. It is this part 
of the trade-off relation, beyond the knee in the information- entropy trade-off, that tells us 
about hypersensitivity to perturbation in the system dynamics. Notice that we need M ^> d 
to investigate this region, but we do not need M so large that a random collection of vectors 
would sample generic vectors, which requires at least M ~ 2 d vectors, i.e., what we call 
a dense collection. Our stochastic perturbation need only produce a sparse collection of 
vectors to see evidence of hypersensitivity; we can say that the vectors in such a sparse 
collection are pseudo-random instead of random (ljj . 

Projective Hilbert space can never be partitioned exactly into spheres of uniform radius. 
This has little effect when the spheres are tiny and numerous, but it becomes a problem when 
there are just a few spheres and prompts us to treat the sphere-grouping trade-off relation 
with caution in this situation. In particular, for the case of just two groups, i.e., J m i n = 1, 
a better method for grouping random vectors is to partition projective Hilbert space into 
two equal volumes defined by the closeness to two orthogonal subspaces of dimension d/2. 
The resulting trade-off for 7 m i n = 1 is analyzed in Appendix [B] and summarized below. We 
abandon the trade-off relation (|2.17j) entirely for J min < 1 because a grouping into spheres 
of uniform radius makes no sense when there are fewer than two spheres. 

We are now in a position to introduce a quantitative measure of a map's hypersensitivity 
to perturbation. For this purpose we introduce the quantity 

1 tt ; (2-18) 



Hs — H(I min = 1) 

1/s is the reduction in system entropy purchased by gathering one (optimal) bit of informa- 
tion about the environment. In a Hilbert-space context, s is an indicator of the randomness 
in a collection of vectors. It is a considerably more informative indicator than the entropy. 
For example, the members of an orthonormal basis together achieve maximal entropy, yet 
a grouping of these vectors into two equally sized groups gives s — 1, independent of the 
dimension. In contrast, for vectors distributed randomly across Hilbert space, the sphere- 
grouping trade-off relation (|2.17|) gives 

1 d 
S = logd - H(I miQ = 1) = (1 + In 2) log(l + In 2) - 1 ~ 3 ' 5 ' (2 ' 19) 

While s provides a signature of randomness, it is the change in s with time, as applied to 
the perturbed system vectors, which indicates the degree to which a system is hypersensitive 
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to perturbation. A rapid increase in s over time has been proposed as a criterion of chaos 
for both classical and quantum systems . In particular, if s increases exponentially with 
time, we say that the system exhibits exponential hypersensitivity to perturbation |g, |8(. 

A detailed analysis p of stochastic perturbations of classically chaotic maps described by 
a symbolic dynamics shows that for such systems, s is indeed a measure of the phase-space 
stretching and folding characteristic of chaotic dynamics. Specifically, s grows as 2 Kt , where 
K is the Kolmogorov-Sinai entropy of the dynamics showing that these systems do 
display exponential hypersensitivity to perturbation and that exponential hypersensitivity 
is equivalent to the standard characterization of classical chaos via the the Kolmogorov- 
Sinai entropy, which in turn is equivalent to characterization in terms of Lyapunov exponents. 

For quantum systems, s is a measure of how fast and how fully the state of a perturbed 
system explores the system Hilbert space. An exponential increase in s indicates both that 
the number of dimensions, d, explored by the perturbed vectors grows exponentially and that 
the vectors populate the explored dimensions randomly. Thus s provides a direct dynamical 
characterization of quantum chaotic dynamics, a characterization that is analogous to the 
characterization of classical chaos in terms of sensitivity to initial conditions. The reason that 
hypersensitivity to perturbation goes beyond the Zurek-Paz chaos criterion of linear entropy 
increase under stochastic perturbations jH 0| is clear: a linear entropy increase indicates 
that the perturbed vectors explore an exponentially increasing number of dimensions, but 
is silent on whether those dimensions are explored randomly. 

A related parameter for characterizing hypersensitivity is the slope of the information- 
entropy trade-off, \dI m i n /dH\, evaluated at J min = (i.e., H = H s ) or perhaps at J min = 1. 
Both the classical analysis in jfjj and the analysis of Appendix [X] prompt us to shy away 
from using the slope evaluated at I m { n = 0, since there are uncertainties about the behavior 
of the slope for very small values of J m i n . Moreover, the slope evaluated at J m i n = 1 seems 
to have no advantages over the parameter s . Thus, in this paper, we calculate numerically 
information-entropy trade-offs for the perturbed quantum baker's maps, and from these we 
determine the time evolution of the hypersensitivity parameter s , preferring it to the more 
problematic use of the slope. 

Having settled on s as our signature of hypersensitivity, we can formulate a better 
information-entropy trade-off for random vectors when 7 min = 1, i.e., for the case of two 
groups. An optimal way of grouping a sufficiently dense collection of random vectors, an- 
alyzed in Appendix [0 is then the following: choose two orthogonal subspaces, each of 
dimension d/2, and partition projective Hilbert space into two equal volumes defined by the 
distance in Hilbert-space angle to these subspaces. The entropy of each partition is 

H = ~ (-A+ log A+ - A_ log A_) « log d ^— (2.20) 



2 V t o -r - / - ndln2 

[cf. Eqs. (EH and (fB23j) with n = d/2], where 

A± = -Jl ± „ Jr/ d L„„ I ~ ~ I 1 + w — I (2.21! 



d V 2 d [{d/2) 



'12 




[cf. Eqs. (fB20jl and (jEZB with n = d/2]. The approximate expressions on the right hold for 
large d and give 

s = - \ — - = 7rdhi2 w 2.2d . (2.22) 

log a — H 
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The coefficient 2.2, smaller than the 3.5 of Eq. ([2. 19)) . indicates that this is a better way to 
partition random vectors into two groups. This value of s represents an approximate upper 
bound for any collection of vectors in Hilbert space. 

Another scenario that is important for the current study occurs when the perturbed 
vectors are restricted to product states of N qubits. Random product vectors can be grouped 
into the two groups corresponding to J m i n = 1 by partitioning the projective Hilbert space 
of one of the qubits into two equal volumes, just as above. The entropy of each partition for 
this qubit is 2 — (3/4) log 3 [Eq. ()2.20|) with d — 2]. Thus for all the qubits, the entropy of 
each partition is 

H = iV-^log3 + l , (2.23) 

which gives 

1 4 

s = = ss 5.3 , (2.24) 

N-H 31og3-4 ' v ' 

independent of D = 2 N . This value represents a rather restrictive approximate upper bound 
on s for product vectors. 

Suppose that for random product vectors, we partition the projective Hilbert spaces of 
j constituent qubits into two equal volumes, thus using j = 7 min bits of information to 
purchase a reduction of the entropy to 

#(/min) = N - Q log 3 - l) J min « N - W5-3 , (2.25) 

for iV > /min > 0. This information-entropy trade-off, which, unlike Eq. (j2.17j) . is linear 
near the maximal entropy, is plotted as the lower dotted line in Fig. |H] (Sec. |VJ). It shows 
that nonentangling quantum maps are not hypersensitive to perturbation. 



B. Other perturbation-based criteria for quantum chaos 

Fidelity decay as a criterion for quantum chaos was introduced by Peres 0,0] (see also [13, 
0,0] and references therein). One compares the unitary evolution of an initial state |^o) 
under the action of a quantum map B with the evolution of the same initial state under 
the action of a modified map, B' — t/B, where the unitary map U is close to the identity 
operator. According to this criterion, a quantum map is chaotic if the fidelity, 

F{t) = \(^\{m^) t B t \^)\\ (2.26) 

decreases exponentially with the number of iterations at short times. In contrast to the 
criterion of hypersensitivity to perturbation, where the effects of a stochastic perturbation 
are analyzed, fidelity decay focuses on just two perturbation histories, corresponding to the 
unperturbed evolution and to a modified evolution where the same perturbation operator 
U is applied at each time step. 

Linear entropy increase as a chaos criterion was introduced by Zurek and Paz 0, 0] . Ac- 
cording to this criterion, a quantum map is chaotic if the entropy ()2.2|) of the reduced system 
density operator (j2.1j) increases linearly with the number of iterations at short times. As we 
have already discussed, a linear entropy increase is essential for exponential hypersensitivity 
to perturbation, but it is not the whole story. 
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III. QUANTUM BAKER'S MAPS 

The classical baker's map is a standard example of chaotic dynamics [j| . It is a symplectic 
map of the unit square onto itself, defined through the equations 

q n+l = 2q n - [2q n \ , (3.1) 
Pn+i = (Pn+ [2q n \)/2 , (3.2) 

where q,p G [0,1), [x\ is the integer part of x, and n denotes the nth iteration of the 
map. Geometrically, the map stretches the unit square by a factor of two in the q direction, 
squeezes by a factor of a half in the p direction, and then stacks the right half onto the left. 

Interest in the baker's map stems from its straightforward symbolic-dynamical char- 
acterization in terms of a Bernoulli shift on binary sequences. If each point of the unit 
square is identified through its binary representation, q = 0-siSz... = YHk=i s k2~ k and 
p = 0-sqS-i . . . = Y^kLo s -fc2 -fc_1 (si G {0, 1}), with a bi-infinite symbolic string 

S = ... S_2S-i3q « SiS 2 53 ■ • • , (3-3) 

V 1 

then the action of the baker's map is to shift the position of the dot by one digit to the 
right, 

(3.4) 

It seems natural to consider a quantum version of the baker's map for the investigation 
of quantum chaos. There is, however, no unique procedure for quantizing a classical map: 
different quantum maps can lead to the same classical baker's transformation. 

To construct a quantum baker's map, we work in a D-dimensional Hilbert space, 7~Cd, 
spanned by either the position states \qj), with eigenvalues qj = (j + 1/2) /D, or the mo- 
mentum states \pk), with eigenvalues pu = (k + l/2)/D (j, k = 0, . . . , D — 1). The constants 
of 1/2 determine the type of periodicity assumed for the position and momentum states, 
in this case, \qj+D) = — \qj), \pu+d) = — \Pk), and thus identify TCd with a toroidal phase 
space with antiperiodic boundary conditions. The vectors of each basis are orthonormal, 
(ijllk) = (Pj\Pk) = 8jk, and the two bases are related via the discrete Fourier transform Fd, 

(qj\F D \q k ) = (qj\Pk) = -L e l ^' h . (3.5) 

V L> 



For consistency of units, we must have 2nhD = 1. 

The first work on a quantum baker's map was done by Balazs and Voros Assuming an 
even- dimensional Hilbert space with periodic boundary conditions, they defined a quantum 
baker's map in terms of a unitary operator B that executes a single iteration of the map. 
Saraceno |2l| later improved certain symmetry characteristics of this quantum baker's map 
by using antiperiodic boundary conditions as described above. To define the Balazs- Voros- 
Saraceno unitary operator in our notation, imagine that the even- dimensional Hilbert space 
is a tensor product of a qubit space and the space of a (D/2)- dimensional system. Writing 
j = x(D/2) + j', x G {0,1}, we can write the position eigenstates as \qj) = \x) <g> \f), 
where the states \x) make up the standard basis for the qubit, and the states \f) are a 
basis for the (D/2)-dimensional system. The state of the qubit thus determines whether the 
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position eigenstate lies in the left or right half of the unit square. The Balazs-Voros-Saraceno 
quantum baker's map is defined by 

B = Fd ° <E> F D ) 2 ) , (3.6) 

where 1.2 is the unit operator for the qubit, and -Fd/2 is the discrete Fourier transform on 
the (D/2)-dimensional system. The unitary B does separate inverse Fourier transforms on 
the left and right halves of the unit square, followed by a full Fourier transform. 

For dimensions D = 2 , an entire class of quantum baker's maps can be defined in 
analogy with the symbolic dynamics for the classical baker's map [101. In this case, we can 



model our Hilbert space as the tensor-product space of iV qubits, and the position states 
can be defined as product states for the qubits in the standard basis, i.e., 

\qj) = \xi) ® \x 2 ) <g) • • • <g> \x N ) , (3.7) 

where j has the binary expansion 

N 

j =x 1 ...x N -0 = Y,xi2 N ~ l (3-8) 
1=1 

and qj = (j + 1/2)/ D — 0-x± . . . xjfl. 

To make the connection with the symbolic dynamics for the classical baker's map, we 
proceed as follows. The bi-infinite strings (J3.3|) that specify points in the unit square are 
replaced by sets of orthogonal quantum states created through the use of a partial Fourier 
transform 

G n = i 2 n <g> F 2 n-u , n = 0,...,N, (3.9) 

where I2™ is the unit operator on the first n qubits and F 2 n-n is the Fourier transform on the 
remaining qubits. The partial Fourier transform thus transforms the N — n least significant 
qubits of a position state, 

Gn ® • • ■ <E> \x n ) <8> K) <g) • • • <8> \a N _ n ) 

= \xt) (8) • • • <g> \x n ) ® ~X= E kn+i) ® ■ ■ ■ ® \x N )e 2max / 2N ~ n , (3.10) 

V 2 x n +l,—,XN 

where a and x are defined through the binary representations a — a\ . . . a^- n - \ and x = 
x n+ i . . . xat-1. In the limiting cases, we have Go = Fd and Gn = The analogy to the 
classical case is made clear by introducing the following notation for the partially transformed 
states: 

I a N _ n . . . ai • xx . . . x n ) = G n \xx) <8> ■ ■ ■ <S> \x n ) <8> \a x ) ® ■ ■ ■ ® \a N _ n ) . (3-11) 

For each value of n, these states form an orthonormal basis and are localized in both position 
and momentum. The state |a7v- n • • • cbi • %i ■ ■ ■ x n ) is strictly localized in a position region 
of width l/2 n centered at 0-xx-- -x n l and is roughly localized in a momentum region of 
width l/2 N ~ n centered at 0-ax . . . ajv-n.1- I n the notation of Eq. ()3.3|) . it is localized at the 
phase-space point la^-n ■ ■ ■ ax • xx . . . x n l. Notice that \a^ ■ ■ ■ = Go |ai) ® ■ ■ ■ (g> |ojv) is 
a momentum eigenstate and that |« Xx ■ ■ ■ xn) = Gn \xx) ® • • - ® \xn) = i\xx) ® • • • ® \xn) 
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is a position eigenstate, the factor of i being a consequence of the antiperiodic boundary 
conditions. 

Using this notation, a quantum baker's map on N qubits is defined for each value of 
n = 1, . . . , N by the single- iteration unitary operator [l0| 

Bn,ti = G n -i ° S n ° G n 1 

= Y I a N-n ■ ■ -a>iXi • x 2 • • .x n )(a N _ n . . . a x • x x x 2 ■ ■ .x n \ , (3.12) 

Xl,...,X„ Ol,— ,OJV-n 

where the shift operator S n acts only on the first n qubits, i.e., S n \xi) <8> \x 2 ) <8> ■ ■ ■ <8> \x n ) ® 
|x n+ i) ® • • • ® \xn) = \x 2 ) <S> ■ • • ® |x n ) ® ® |x n+ i) ® • • • <S> \xn). Notice that since S n 
commutes with G~ l , we can put B^ n in the form 

B N>n = i 2 «-i ® (F 2 ^-„ + i o (i 2 <gi #^_ n )) o »S n . (3.13) 

Since Si is the unit operator, it is clear that 5jv,i is the Balazs-Voros-Saraceno quantum 
baker's map (J3.6|) . It is worth mentioning here that Ermann and Saraceno 2^1 have recently 



proposed and investigated an even larger family of quantum baker's maps, which includes 
all of the above quantizations as members. For the purposes of this article, however, we 
need only consider B^ yn . 
We can also write 

Bn,u = 2-2™- 1 ® Sjv-n+l,l ° S n , (3-14) 

which shows that the action of i?jv,n is a shift of the n leftmost qubits followed by an 
application of the Balazs-Voros-Saraceno baker's map to the N — n + 1 rightmost qubits. 
At each iteration, the shift map S n does two things: it shifts the nth qubit, the most 
significant qubit in position that was subject to the previous application of Sjv-n+i,i, out 
of the region subject to the next application of B N _ n+1>1 , and it shifts the most significant 
qubit in position (first qubit) into the region of subsequent application of B N _ n+ i i. 

The quantum baker's map B^^ takes a state localized at lazv-n • • • Qi • %i ■ ■ - x nX to a 
state localized at laAr_ n . . . a\X\ *x 2 . . . x n l. The decrease in the number of position bits and 
increase in momentum bits enforces a stretching and squeezing of phase space in a manner 
resembling the classical baker's map. In Fig. HJa), (b), (c), and (d), we plot the Husimi 
function (defined as in j23|) for the partially Fourier transformed states (j3.11|) when N = 3, 
and n = 3, 2, 0, and 1, respectively. The quantum baker's map is a one-to-one mapping of 
one basis to another, as shown in the figure. 

One useful representation of our quantum baker's maps, introduced in [T^j, starts from 



using standard techniques 24] to write the partially transformed states ()3.10|) as product 
states: 



a N _ n . . . ai • xi . . . x n ) 



n 

3 7ri(0-oi...ajv-n 1 ) 




These input states are mapped by B Nn to output states 
| a N - n ■ ■ ■ aiXi •x 2 ...x n ) — e 



( ® -J=(|0) + e 2 ^ ( °- ajv - fc + 1 - ajv -™ 1) |l))^ . (3.15) 
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(a) 



(b) 




(c) 



(d) 



FIG. 1: Husimi function for each partially Fourier transformed state 1 If) when N = 3: (a) n = 3, 
(b) n = 2, (c) n = 0, and (d) n = 1. The action of the quantum baker's map is to map the 
eight states in (a) to the eight states in (b), as shown by the numbers labeling the states. Similarly, 
and -E^i map one set of partially Fourier transformed states to another, as indicated by the 
arrows. The map B% i is the Balazs-Voros-Saraceno quantum baker's map. 



/ N i \ i 

<g>( (g) (|o) + e 2 " (0a ^- fe + 1 - ajv -" 1) |l)) ® (| ) + e 2 "(°- xiai - a ^-" 1 )|l)) . (3.16) 

These forms show that the quantum baker's map Bn,u shifts the states of all the qubits to 
the left, except the state of the leftmost qubit. The state of the leftmost qubit can 
be thought of as being shifted to the rightmost qubit, where it suffers a controlled phase 
change that is determined by the state parameters a\, . . . , a^-n of the original "momentum 
qubits." The quantum baker's map can thus be written as a shift map on a finite string 
of qubits, followed by a controlled phase change on the least significant qubit in position. 
In |25( this shift representation was developed into a useful tool. Using an approach based 
on coarse graining in this representation, the classical limit of the quantum baker's maps 
was investigated. 

The classical limit for the above quantum baker's maps was also investigated in |23| . 
using an analysis based on the limiting behavior of the coherent-state propagator of B^^. 
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When D = 2 — > oo, the total number of qubits N necessarily becomes infinite, but one has 
a considerable choice in how to take this limit. For example, we could use only one position 
bit, thus fixing n = 1, and let the number of momentum bits N — 1 become large. This 
is the limiting case of the Balazs-Voros-Saraceno quantization. There is, however, a wide 
variety of other scenarios to consider, e.g., n = N/2 or n = 2N/3 — 1 as iV — > oo. In Q 
it was shown that provided the number of momentum bits N — n approaches infinity, the 
correct classical behavior is recovered in the limit. If the number of momentum bits remains 
constant, i.e., n = N — k (k constant) as N —>■ oo, a stochastic variant of the classical baker's 
map is found. In the special case n = N this variant takes the form 

s = . . . s_ 2 S-iSq • sis 2 s 3 . . . — > s' = . . . s_ 2 s~is r • s 2 s 3 . . . (3.17) 

when written in the symbolic-dynamical language of Eq. (J3.3)) . The bit r takes the value si 
with probability cos 2 [7r/2(0 • s s_is_ 2 • • • — 1/2)1 (and 1 — si otherwise). These results are 
consistent with those obtained previously [25J. 

The extremal map has other interesting properties. All finite-dimensional unitary oper- 
ators are quasi-periodic; the quantum baker's map Bn,n, however, is strictly periodic, 

(B NtN ) = 1 , (3.18) 

as we show below. All its eigenvalues, therefore, are 4N-th roots of unity, i.e., of the form 
e mk / 2N ; a nd hence, there are degeneracies when N > A. This represents a strong deviation 
from the predictions of random matrix theory [26|. The eigenstates of the extremal map 
were recently studied by Anantharaman and Nonnenmacher [27}, where B^,n (with periodic 
rather than antiperiodic boundary conditions) was called the "Walsh-quantized" baker's 
map. The above degeneracy in the eigenvalues allows constructions of eigenstates that 
remain partially localized in the semiclassical limit, which means that "quantum unique 
ergodicity" |28[ fails for this quantization. 

The periodicity of the extremal map ()3.18|) can be easily shown after noting that Bn,n = 

—iGn-i Sn = — «(i-2 JV - 1 <8> F2) o Sn] i.e., Bn,n is a shift followed by application of the 
unitary 

C = -iF 2 = -L( e -**/ 4 (|o}(0| + |1)(1|) + e m/4 (|0)(l| + |1)(0|)) = e --/4 e (W4)^ ; (3 19 ) 

which is a rotation by 90° about the x axis, to the least significant position qubit. On 
product states, the action of B N ^ N can be written explicitly as 

B n ,n\^i) ® ^2) ® • • • ® iV'iv) = \ifa) ® • ■ • ® |V>iv) ® C'lV'i) • (3-20) 

Since C 4 = 1, we get the property ()3.18|) . One can also see that B N)N cannot entangle initial 
product states. 

When n < N, the action of the quantum baker's map is similar to Eq. ()3.2()j) . but 
with a crucial difference. After the qubit string is cycled, instead of applying a unitary to 
the rightmost qubit, a joint unitary is applied to all of the N — n + 1 rightmost qubits. 
As discussed above, this joint unitary can be realized as controlled phase change of the 
rightmost qubit, where the control is by the state parameters a 1; . . . ,aAr_ n of the original 
momentum qubits. This controlled phase change means that initial product states become 
entangled. 
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Entanglement production under iterations of the quantum baker's maps was the subject 
of a recent paper j^. Since the entangling controlled-phase change involves an increasing 
number of qubits as n decreases from n = N to n = 1 (the Balazs-Voros-Saraceno map), one 
might expect that the entanglement increases as n ranges from N to 1. What was found, 
however, is that provided n is not too close to N, all the maps are efficient entanglement 
generators, but the greatest entanglement is produced when n is roughly midway between N 
and 1. Starting with a uniform distribution of initial product states, the mean entanglement 
"quantum-baked" into the distribution was found to saturate at a level near to that expected 
in random states. The small deviations from the entanglement of random states might be 
due to hidden symmetries in the quantum baker's maps |30j | . 

Lastly, we mention another difference between the extremal quantum baker's map -Bjv.at 
and other members of the baker's map family of quantizations. Ermann, Paz, and Sara- 
ceno |3jj have found that when a system with the dynamics of a quantum baker's map is 
cast in the role of an environment acting on another quantum system, the extremal quan- 
tum baker's map Bn,n is l ess effective at inducing decoherence than other members of the 
family. In particular, they showed that while the entropy production rates of the different 
quantum baker's maps are indistinguishable on a short time scale, which scales linearly with 
N, B N> x saturates much sooner than the other maps, thus displaying the behavior expected 
for regular systems. 

In view of the above described anomalous behavior experienced by the extremal map, 
Bn,n, our curiosity now invites an investigation into the various currently prevailing 
perturbation-based tests for quantum chaos, as applied to our class of quantum baker's 
maps. We start, however, by investigating the simplest example, the extremal map Bn,n 
itself. 



IV. CHAOS IN THE EXTREMAL QUANTUM BAKER'S MAP? 
A. The extremal map and perturbations 

In Sec. 11111 we considered different quantizations of the baker's map as unitary transfor- 
mations, B NtU {n — 1, . . . , N), on a set of N qubits. When written in the form (|3.14|) . each of 
these transformations consists of two steps: a cyclic shift, S n , in which the n leftmost qubits 
are shifted without otherwise being altered and a unitary transformation on the rightmost 
N — n + 1 qubits. For the extremal quantum baker's map, n = N, which we consider in 
this section, this second transformation is the gate C of Eq. ()3.19|) . which acts only on the 
single rightmost qubit as in Eq. ()3.20|) and rotates it by 90° about the x axis. 

We first examine the behavior of B^^ under perturbations after each iteration. In an 
effort not to affect qualitatively the dynamics of the map itself, we choose our perturbations 
to be correlated across the smallest possible distances in phase space. One choice might 
then be to perturb only the single rightmost, least significant qubit in the position basis. 
Indeed, such a choice leads to the smallest changes in position. As a consequence of the 
uncertainty principle, however, perturbing the least significant qubit in the position basis 
causes correlated changes across the greatest distances in momentum. In opposition to our 
classical intuition, no single qubit can be thought of as being "more significant" than another 
in an overall phase-space sense. Our particular choice of qubit upon which to perturb does 
not affect the phase-space area of correlated changes made to a state; however, perturbations 
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affecting the middle qubit(s) give rise to correlated changes across the smallest phase-space 
distances. In the present analytical study, it is simplest to take the rightmost qubits as being 
the least significant. These considerations are revisited later in our numerical investigations, 
where we instead choose to perturb the middle qubit(s). 

Suppose that with each step we perturb the m rightmost qubits, where m C A^, by- 
applying an m-qubit unitary transformation chosen at random. For the moment, we 
assume that any transformation is allowed, but the arguments still work even if only a finite 
set of transformations is allowed. 

For simplicity, we assume that the system is initially in the tensor-product state \0)® N . 
Suppose that m = 1 and that the perturbation affects only the single rightmost qubit. Then 
the perturbation operators are all of the form fi^- 1 <8> U% , and the state after a single step 
becomes 

| )®* _> lO)® 7 ^ 1 ® (u£?C\0)) , (4.1) 
while the next step transforms it to 

lO)^- 1 ® (tffftflO)) - |0>^ 2 ® (U^C\0)) ® (U£C\0)) , (4.2) 

and so forth. It is clear that the above dynamics does not explore the entire Hilbert space, 
since the state remains a tensor product as long as the perturbation is restricted to a single 
qubit. 

The situation changes if we let m = 2. Since two-qubit gates between nearest neighbors 
are sufficient for universal quantum computation, any state can be produced by the shift 
map plus two-qubit perturbations. It does not follow, however, that all states can be reached 
quickly; in general, the number of gates needed to reach a generic state of N qubits increases 
exponentially with N, which implies that many "rounds" (complete sets of N steps) are 
needed to reach most states. On the other hand, as we saw in Sec. Ill A 21 the perturbation 
need not sample generic vectors to elicit evidence for hypersensitivity, so considerations of 
universality in quantum computation and the time needed to sample generic states provide 
little information about hypersensitivity. 

B. Signatures of chaos for Bn,n 

In this section we show that the extremal quantum baker's map is chaotic according to 
two popular signatures of quantum chaos: it displays a linear entropy increase when coupled 
to an environment, and the fidelity between two vectors evolving according to the original 
map and a slightly changed version of the map decreases exponentially. It does not, however, 
display hypersensitivity to perturbation. 

1. Fidelity decay 

Define a modified baker's map by 

Bn,n\ x i) ® • • • ® \ x n) = ® • • • ® \x N ) ® U (1) C\xi) , (4.3) 

where is a single-qubit unitary map satisfying < | (OlC^U^ClO) \ < 1. We can define 
A = -21n|(0|(7tf>( 1 )C'|0)| > 0. Letting |^(t)) = {B NtN f\ip ) and \^{t)) = 
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where \if) ) = \0)® N , we see that the fidelity decreases exponentially with the number of 
iterations" 

F(t) = \W(t)mt))\ 2 = e- Xt , t = 0,...,N. (4.4) 

For greater numbers of iterations, the simple exponential decay is modified as the qubits 
experience more than one application of C or U^C. 

2. Linear increase of entropy 

Let the environment be a collection of qubits in the maximally mixed state. After each 
iteration of the map, the register interacts with a fresh environment qubit. The interaction 
is given by a controlled a x operation, with the environment qubit acting as control and 
the target being the rightmost system qubit. In the notation of Sec. Ill A 1\ this binary 
perturbation amounts to an application of one of Uq = 1 2 jv or U\ = fi^jv-i ® &x (chosen 
with equal probability) at each time step (this stochastic perturbation is the m — 1 model 
of Sec. II V A| with the single-qubit perturbation unitaries restricted to the identity and a x ). 
After tracing out the environment, one iteration of the perturbed map is described by the 
quantum operation 

B{p) = ]^B N)N pB ] N N + ^ (i 2 w-i ® Si) B N , N pB ] NN (i 2 w-i ® a x .) . (4.5) 

For the initial state p = (\0)(0\)® N , and denoting by £>* the t-th iterate of B, we have 

B*(A)) = (|0><0|) 8 ^-*)® (i 2 /2) W , (4.6) 

since a x commutes with C. The entropy of <B*(po) is t bits. The entropy thus increases at 
a rate of 1 bit per iteration until it saturates at N bits after N steps. Under single-qubit 
perturbations, it is clear that the perturbed vectors explore an exponentially increasing 
number of Hilbert-space dimensions, but it is equally clear that they do not explore these 
dimensions randomly. 

3. Hypersensitivity to perturbation 

From the discussions in Sees. Ill A~2l and HV Al it can be seen that Bn,n is not hypersensi- 
tive to perturbations that affect only the single rightmost qubit. For an extreme example of 
this, consider the binary perturbation in Sec. IIV B 21 immediately above. After t < N map 
iterations, there are 2* perturbation histories, which, with the initial state |^o) — |0)® , 
correspond to orthogonal system vectors: 

\^) = \0)® {N - t) C\k 1 )®---®C\k t ) , (4.7) 

where k{ G {0, 1}, in the notation of Sec. Ill A II A measurement on the environment that 
groups these vectors according to the values of ki,...,kj, with < j < t, reduces the 
average system entropy from H s = t to H(I min — j) = t — I min bits. When t = N, the 
perturbed system vectors make up an orthonormal basis, and for t > N, the perturbation 
produces 2 l ~ N copies of an orthonormal basis. Thus, for t > N, the information-entropy 
trade-off relation is H(I min ) = N — J min . For all t, our hypersensitivity parameter takes the 
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value s = 1. In this extreme example, each bit of information purchases a bit of entropy 
reduction, as is always true when the perturbed vectors are drawn from an orthonormal 
basis with each vector in the basis having the same overall probability. 

In general, stochastic perturbations that affect only a single qubit of the extremal quan- 
tum baker's map are expected to produce an information-entropy trade-off that is linear 
(to good approximation) near the maximal entropy. Although the hypersensitivity param- 
eter generally varies with both the choice of perturbation and number of map iterations, 
its magnitude should not exceed 5.3, the bound on s for product states [Eq. (|2.24j) ]. The 
extremal map -Bjv,7v, therefore, does not exhibit exponential hypersensitivity to perturba- 
tion under single-qubit perturbations. By contrast, we have seen above that it does exhibit 
linear growth of entropy and exponential decay of fidelity. Hypersensitivity to perturbation 
is evidently a finer sieve than the other two perturbation-based criteria. 

The reason the perturbed extremal map does not explore Hilbert space efficiently is 
that the map itself produces no entanglement. In contrast, the nontrivial quantizations of 
the baker's map are efficient entanglement generators |29j . producing entanglement that 
saturates after several iterations at a level close to that expected in random states. For 
these nontrivial quantizations, even a single-qubit perturbation, together with the entangling 
transformation of the unperturbed map, generically gives rise to a universal set of unitary 
gates, so in time the system can approach any state in the Hilbert space. Although the 
speed at which this happens remains unknown, our numerical results for hypersensitivity 
to perturbation, presented in the next section, suggest that if n is not too close to N, the 
perturbed nontrivial quantizations do efficiently explore all of Hilbert space. 

Both the simple analytical argument above and the numerical results in the next section 
are for single-qubit perturbations. A systematic study of hypersensitivity to perturbations 
acting on two or more qubits is beyond our current numerical capabilities. In the remainder 
of this section, we present an analytical argument that suggests that, for a small number of 
time steps and for maps close to the extremal map Bn^Ni the information-entropy trade-off 
is linearly bounded even for entangling perturbations acting on two qubits. 

We choose a perturbation that affects the two rightmost qubits, i.e., k = 2. Given some 
reasonable assumptions about the stochastic perturbation, if we average over all perturba- 
tions, the state after t steps is approximately equal to 

pW^(|0)(0|)^- t - 1 ®(i 2 /2)^ +1 . (4.8) 

This state has von Neumann entropy H$ = t + 1. We would like to acquire some information 
J m i n about the perturbations which enables us to reduce this entropy by a small amount 
AH = H s - H. 

We now show that the ratio J m in/ AH is bounded above by a quantity that is independent 
of t for all t < N. Suppose that after t steps, our system is in state (|4.Bp . Now let us apply 
B n ,n, but not the perturbation. If we trace out all but the two least significant qubits, these 
two qubits are in the state 

p® = i 2 /2®6\0)(0\& , (4.9) 

which has one bit of entropy. The perturbation affects only these two bits, so the state of 
the other N — 2 qubits is irrelevant to the entropy increase. Now we apply the perturbation 
and get 

P (2) -E^P (2) ^ = (W2)^ , (4.10) 
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where £ labels which perturbation is performed, p^ is the probability of that perturbation, 
and is the corresponding two-qubit unitary transformation. The entropy of the new state 
is two bits, giving an entropy increase of one bit. 

Clearly, we can reduce the entropy by one bit if we can determine which was actually 
performed. If the perturbations are drawn from a discrete set, the number of bits needed to 
determine this is given by the entropy of the distribution i.e., 

If £ is continuous, then fully determining would require an infinite amount of information. 
Since the space of two-qubit operators is not very large, however, it doesn't take that much 
information to know to a good approximation; e.g., we could achieve an entropy reduc- 
tion of nearly a bit at a cost of approximately 45 bits by knowing each of the 15 relevant 
parameters of an arbitrary two-qubit unitary with three bits of precision. 

This procedure, while not necessarily optimal, places a rather low bound on the ratio 
Jmin/Aff, a bound independent both of the number of iterations, t, and the number of 
qubits, N. This argument changes little if we use B NjN _i instead of B NjN , or B NjN _ k for k 
small compared to N. Nor does it change much if the perturbation affects k bits, so long as k 
is small compared to N. If the perturbation affects many bits, however, or if a quantization 
Bn,n-u is used for large k, the upper bound on I min /AH becomes so large that it gives little 
restriction. 

The above provides some evidence for the conjecture that maps close to the extremal 
map, B N>Nl do not exhibit exponential hypersensitivity to entangling perturbations. Since 
these results are valid only as long as t < N, i.e., as long as the number of iterations does not 
exceed the number of qubits, this evidence must be regarded as suggestive, but inconclusive. 



V. NUMERICAL RESULTS 

We now investigate numerically the entire class of quantum baker's maps, B^, n ( n — 
1, . . . ,N), in the context of the three perturbation-based criteria for quantum chaos. As 
remarked in Sec. IIV A"l perturbations affecting the middle qubit(s) cause correlated changes 
to the state in phase space across the smallest possible distances. Up until now we have 
applied all perturbations to the rightmost, least significant qubit in position. Since the 
application of C' as a perturbation to this qubit would undo the dynamics in momentum in 
the case of the extremal quantum baker's map, B N ^ N , one might judge this perturbation to 
be atypical, upsetting the crucial momentum dynamics of the map. To avoid this, we choose 
henceforth the total number of qubits N to be odd, and we perturb the middle qubit. 

The perturbation we choose for this qubit is a simple binary perturbation, a rotation by 
angle ±27ra about the y axis, 

U k (a) = i 2( iv-i)/2 <g> e li(4) '^ <g> i 2( iv-i)/2 . (5.1) 

The perturbation is conditioned on the binary environmental states |A;}#, k — or 1. To 
be precise, after each iteration of the map, the system couples to its environment through a 
joint conditional evolution with end result 

Ptotai = \ [(UopUS) ® |0) B <0| + (UipUl) ® |1M1|] • (5.2) 
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FIG. 2: The rate of increase in entropy is initially linear for all (a) five-qubit and (b) eleven-qubit 
quantizations of the baker's map. The rate of decrease in fidelity is initially exponential for all (c) 
five-qubit and (d) eleven-qubit quantizations. 



To avoid inapt comparisons, we use a = 0.2 (rotation angle OAir) and initial system state 
IV'o) — |O) 0Ar throughout this section. 

We can cope with the unwanted perturbation in many different ways. One possibility 
is to accept an increase in entropy and average over all perturbation histories by tracing 
out the environment. The system entropy, H$ = — tr(plogp), then increases at an initially 
constant linear rate for all quantum baker's maps. This is shown for the quantizations using 
N = 5 qubits in Fig. El^a) and N = 11 qubits in Fig. |2^b). The rate of entropy production 
for the different quantizations is nearly the same for the first t = (N — l)/2 iterations. 
In contrast to the other quantizations, the entropy produced by the extremal map, Bn,n, 
remains constant for times beyond t = (N — l)/2 before resuming its climb towards the 
maximal entropy of 5 or 11 bits. Thus, although there is a quantitative change in entropy 
production at later times, the different baker's maps behave qualitatively the same. These 
numerical results support the simple analysis of Sec. IIV B 21 

Alternatively, if the above entropy production proves unacceptable, we could instead 
perform a measurement on the environment at each time step to record which perturbation 
actually occurs. Consider the perturbation history of all l's. The fidelity decay between two 
initially equal quantum states that evolve either according to this extremal perturbation 
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history or the unperturbed map, 

(5.3) 

might be used as an indicator of the underlying dynamics of the map. The rate of fidelity 
decay for all quantum baker's maps is initially exponential. This is shown in Figs. |2f c) and 
I2td) for the five-qubit and eleven-qubit quantizations. Although the fidelity corresponding to 
the extremal map, Bn,n, stalls at approximately t = (N + l)/2 iterations, all quantizations 
are found to exhibit decay rates which are initially exponential. Again, these numerical 
results support the simple analysis of Sec. IIVB 11 

The iteration at which the entropies and fidelity decays first become appreciably different 
for the various quantizations remains at t — (N + l)/2 as N increases, and thus our conclu- 
sions become stronger in the limit of large N. To keep our analysis strictly in the quantum 
regime, however, we focus on the five-qubit quantizations for the remainder of this section. 

To investigate hypersensitivity to perturbation, we first consider a particularly intuitive 
algorithm for grouping vectors, which is based on finding structure produced by the temporal 
order of the perturbations. Each grouping corresponds to measuring, after a fixed number of 
iterations t, the environment states — and, hence, the applied perturbation — at I < t times. 
The 2* perturbation histories — and their final states — are thus grouped into 2 l sets, each 
containing 2 f ~ l states. It takes I = I bits to specify a group. 

As an example of this procedure, suppose we have t = 4 iterations and we choose to 
measure the first and last states of the environment. Thus I = 2, and all histories are 
grouped into 2 l = 4 sets of 2 t ~ l = 4 binary strings in the form 0**0, 1**0, 0**1 and 
1**1, where * denotes an arbitrary entry. Defining 

\kxk 2 ■ • ■ kt) = U kl B Ntn U k2 B Ntn ■ ■ ■ U kt B Njn \ip ) , (5.4) 

the final state of the system, conditioned on measurement results % and j for the first and 
last environment qubits, is 

Pij = 2 l ~ t Yl \ik2k 3 j)(ik 2 k 3 j\ . (5.5) 

k2,k 3 &{0,l} 

Consequently, at the expense of storing I = 2 bits of information, we can, on average, reduce 
the entropy to 

H = -7Ti J2 tr (pij log fa j) . (5.6) 

i,je{o,i} 

The particular two bits stored in this example might not be the optimal choices. There 
are countless other measurements to consider, some of which no doubt lead to lower average 
entropies. For the moment, however, we restrict our measurements to the above type and 
minimize H over the (V\ possible choices for the measurement times. Denoting this minimum 
entropy by H, the minimum information needed to reduce the average system entropy to H is 
then J m i n = I bits. Although there is no guarantee that I is in fact the overall minimum, this 
simple scheme, which we call the temporal grouping algorithm, proved superior to previously 
used schemes (e.g., those discussed in 0) for the maps and perturbations considered here. 

Using this procedure, in Fig. El we plot 7 m i n versus H for all perturbed five-qubit quantum 
baker's maps after t — 16 map iterations (solid lines). The perturbing parameter remains 
at a = 0.2 and initial state at \ipo) = (0)®^. The region of interest regarding the question 
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FIG. 3: The minimum information I m { n needed to reduce the entropy to H after 16 map iterations, 
using the temporal grouping algorithm, for all perturbed five-qubit quantizations of the baker's 
map (n = 1, . . . , 5). The upper and rightmost dotted curves are the approximate sphere-grouping 
trade-off for random vectors in d = 32 dimensions [Eq. (|2.17|) ]. while the lower dotted line is 
the linear trade-off for random product vectors [Eq. ()2.25|) ]. The inset shows the hypersensitivity 
parameter s for each quantization. 

of quantum chaos lies to the right where I m \ n is small. Here we see that, except for the 
quantizations with n close to N = 5, a very large amount of information is required to 
reduce the system entropy by a small amount. This is a distinguishing characteristic of chaos, 
which is absent for the extremal quantization (n = N). When I min is small, the information- 
entropy trade-off is characterized by the hypersensitivity parameters [Eq. ([2.18)1 ] . The inset 
in Fig. El shows this quantity for all five quantizations. Recall that 1/s is the reduction 
in system entropy purchased by gathering one bit of information about the environment. 
Although entropy reduction is affordable for the extremal quantization, one bit buys very 
little when n approaches 1. The dotted lines show our theoretical trade-offs for random 
vectors (upper and rightmost) and random product vectors (lower), given by Eqs. ([2.17)1 
and ()2.25[) . respectively. When n approaches 1, the information-entropy trade-off approaches 
that expected for random vectors, while for n = N, it is bounded by the trade-off for random 
product vectors. 

Using the same grouping algorithm, we plot in Fig. 0] the information-entropy trade-off 
for a growing number of iterations of the Balazs-Voros-Saraceno quantization (n = 1). The 
figure shows J m i n versus H for 5-18 iterations of -85,1, and in the inset, the corresponding 
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FIG. 4: The minimum information I m { n needed to reduce the entropy to H after t map iterations, 
where 5 < t < 18, for the perturbed five-qubit Balazs-Voros-Saraceno quantization of the baker's 
map (n = 1), using the temporal grouping algorithm. The dotted curve is the sphere-grouping 
trade-off for random vectors in 32 dimensions. The inset shows the hypersensitivity parameter s 
at each iteration. 

value of s. To a rough approximation, our hypersensitivity signature s appears to grow 
exponentially with the number of iterations. This map thus exhibits numerical evidence 
of exponential hypersensitivity to perturbation. Notice, however, that the trade-off violates 
the sphere-grouping bound (|2.17|) derived from random states for t ^ 16 (dotted line), which 
means the current method for gathering information about the environment is not optimal, a 
situation we discuss further below. The extremal quantization displays a strikingly different 
behavior. The information-entropy trade-off and the parameter s for 5-18 iterations of B 5 5 
are shown in Fig. |5J In this case the information-entropy trade-off remains approximately 
linear for all levels of iteration, with a very roughly constant s. There is no evidence of 
hypersensitivity to perturbation. 

We now investigate the hypersensitivity parameter in greater detail for the Balazs-Voros- 
Saraceno and extremal quantizations. The graphs of s in the insets of Figs. H] and 03 are 
redrawn, now on a logarithmic scale, as the dashed lines in Figs. Efa) andEfb), respectively. 
The horizontal dotted lines in each of these figures are the upper bounds ()2.22|) and (|2.24|) . 
respectively, corresponding to the values of s for random vectors and random product vectors. 
Notice that in both cases the dashed lines cross these bounds. This indicates that the 
temporal groupings used up until now are not optimal. Indeed, in the case of the Balazs- 
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FIG. 5: The minimum information / m j n needed to reduce the entropy to H after t map iterations, 
where 5 < t < 18, for the perturbed five-qubit extremal quantization of the baker's map (n = N), 
using the temporal grouping algorithm. The dotted line is the trade-off for random product vectors 
[Eq. (|2.25|) ]. The inset shows the hypersensitivity parameter s at each iteration. 

Voros-Saraceno quantization, by considering groupings which correspond to partitions of 
projective Hilbert space into two equal volumes, we find that one bit of information can buy 
larger entropy reductions when t > 16. This method of grouping vectors, however, works 
well only for distributions that are close to random. 

We now consider grouping algorithms that are not constrained by a supposed tempo- 
ral structure of the vector distribution. Although optimal groupings can always be found 
by simply testing every possibility, the size of the search space is doubly exponential in the 
number of map iterations. We thus turn to the theory of combinatorial optimization. Specif- 
ically, a simple genetic algorithm [32| was used to partition vectors into two groups with 
the goal of minimizing the average conditional entropy 77 Although these groups were not 
constrained to be of equal size, the returned solution always corresponded to 7 = 1 ± 0.003, 
and thus we can take H(I min — 1) = H to a very good approximation. The corresponding 
value of s for this method is plotted as the solid line in Fig. El In many cases the genetic algo- 
rithm located precisely the same vector grouping that was found previously by the temporal 
grouping algorithm. The spurious higher values of s , however, are now significantly reduced 
for both quantizations. Although Fig. E(a) remains incomplete due to computational con- 
straints, for the data points calculated, logs has regained its linear approach to the upper 
bound, where it eventually will saturate. The difference between the two quantizations un- 
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FIG. 6: The hypersensitivity parameter s after t map iterations, where 5 < t < 18, for the 
perturbed five-qubit (a) Balazs-Voros-Saraceno (n = 1) and (b) extremal (n = N) quantizations 
of the baker's map. The spurious higher values of s arising from the temporal groupings (dashed 
lines) are significantly reduced using a genetic-algorithm approach (solid lines). The dotted lines 
are the approximate upper bounds on s corresponding to (a) random vectors [Eq. (|2.22[) ] and (b) 
random product vectors [Eq. 1)2.24(1 ], 



der single-qubit perturbations is now difficult to dispute. The criterion of hypersensitivity 
to perturbation thus unmistakeably distinguishes the dynamics of the extremal quantum 
baker's map (n = N) as qualitatively different from the Balazs-Voros-Saraceno quantization 
(n = l). 



VI. CONCLUSION 



This paper addresses the difficult question of how to characterize quantum chaos dy- 
namically in the hard quantum regime, far from the classical limit where signatures of the 
classical sensitivity to initial conditions can be identified in the quantum properties of a 
system. In this hard quantum regime, criteria for quantum chaos rely on studying the 
effects of perturbing the quantum dynamics. In this paper we study the three perturbation- 
based criteria that have been proposed: linear increase of entropy when a system is coupled 
to a perturbing environment; exponential decay of the fidelity between the unperturbed 
dynamics and a modified dynamics; and hypersensitivity to perturbation under stochastic 
perturbations of the dynamics. Hypersensitivity to perturbation is formulated in terms of 
the entropy reduction achieved by acquiring information about the perturbation, which we 
call the information-entropy trade-off. Of these three criteria, hypersensitivity to pertur- 
bation is by far the most difficult to define rigorously and to investigate analytically and 
numerically. 

We apply these three criteria to a set of qubit-based quantizations of the baker's map. 
These quantizations range from a map that has a trivial unentangling shift dynamics to the 
original Balazs-Voros-Saraceno quantized map, which is highly entangling. We find, through 
a combination of analytical arguments and numerical results, that all the quantizations 
exhibit a linear entropy increase and an exponential fidelity decay. In contrast, we show that 
the criterion of hypersensitivity to perturbation distinguishes the entangling quantizations 
from the shift map. In particular, by focusing on the trivial shift map and the Balazs-Voros- 
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Saraceno map, our numerical work on hypersensitivity provides compelling evidence that 
these maps behave quite differently under stochastic perturbations, as revealed by studying 
the information-entropy trade-off for these maps. 

The reason that hypersensitivity to perturbation is different from the other two 
perturbation-based criteria is not hard to identify Linear entropy increase and exponential 
fidelity decay both tell one about how widely a perturbation disperses vectors in Hilbert 
space, but they provide no information about how randomly the perturbed vectors popu- 
late Hilbert space. In contrast, the randomness of the distribution of perturbed vectors is 
precisely what the information-entropy trade-off is sensitive to. 

One way to quantify the trade-off in a single number is provided by the hypersensitivity 
parameter s: s^ 1 is defined in terms of the information-entropy trade-off as the entropy 
reduction purchased by an optimal one bit of information about the perturbation, but s 
can be interpreted as the number of Hilbert-space dimensions explored randomly by the 
perturbed vectors. Thus when s increases exponentially, as our numerical work indicates 
for the Balazs-Voros-Saraceno map, it signals that the perturbed vectors are populating an 
exponentially increasing number of dimensions in a random way. Linear entropy increase 
and exponential fidelity decay do not provide information about this property of chaotic 
quantum dynamics. 

The numerical hypersensitivity results in this paper are obtained for single-qubit per- 
turbations of the five-qubit baker's maps. Investigating more general perturbations would 
involve dealing with more qubits and thus would require considerably greater computa- 
tional resources. To reduce the required computational resources in future investigations of 
hypersensitivity, it would be highly desirable to have an analytical argument or sufficient 
numerical evidence to demonstrate convincingly that the hypersensitivity parameter s is, 
by itself, a reliable signature of hypersensitivity to perturbation. Were this established, 
numerical investigations of hypersensitivity could be reduced from computing the entire 
information-entropy trade-off to calculating the trade-off only for the case of one bit of 
acquired information. 
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APPENDIX A: INFORMATION-ENTROPY TRADE-OFF FOR RANDOM 

VECTORS 

Consider M state vectors distributed randomly in a d- dimensional Hilbert space, where 
we assume that M > d. Given an entropy H < logd, we group the vectors into groups 
that on average have this entropy and then ask how much information is required to specify 
a group. For a given H, we are interested in the grouping that minimizes the required 
information, as in Eq. (|2.9|) . The relation between J m i n and H is the information- entropy 
trade-off. In this Appendix we formulate an approximate information-entropy trade-off for 
random vectors by grouping the vectors into spheres on projective Hilbert space whose radius 
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is given by a Hilbert-space angle 0. This being an approximate trade-off relation, we denote 
the information by J instead of J m i n . 

The sphere-grouping model is based on results, given in , for the volume and entropy 
of a Hilbert-space sphere. The model was formulated in 0] and refined in 11 . 



The number of spheres of radius that can be packed into projective Hilbert space is 
given by Eq. (A18) of Q, 

= = (sin^r^ , (ai) 

where Vd(<t>) is the volume of a sphere of radius an d Vd is the total volume of projective 
Hilbert space in d dimensions. The entropy of a mixture of vectors distributed uniformly 
within a sphere of radius is given by Eqs. (B5)-(B6) of 

H d (<f>) = -A log A - (1 - A ) log(^y) = H 2 (X ) + (1 - A ) log(d - 1) , (A2) 
where H 2 (Xq) is the binary entropy corresponding to the largest eigenvalue 

A = 1 - — sin 2 . (A3) 

If we group the M vectors into groups of radius 0, the number of vectors per group is 

A^(0) = ^y=AT(sin 2 0) d - 1 , (A4) 

provided this number is not less than one. There is a critical angle, 06, at which there is 
only one vector per group, i.e., (sin 2 <p b) dl = 1/Af. For > 06, there are Md{4>) groups, 
each containing approximately AV(0) vectors, but for < 6 , there are M groups, each 
containing just one vector. The information required to specify a group at resolution angle 
is thus 7(0) = log A/" for < 06 and 7(0) = logA/" d (0) for > <p b . There is another critical 
angle, 4> d , at which there are only two groups, i.e., (sin 2 0d) a! ~ 1 = 1/2. For > 0^, we can 
not talk about grouping the vectors into spheres of equal radius, so we remove these angles 
from consideration. Thus we write the information to specify a group as 

JW -\logA' (i (0) = -(rf-l)log(sin 2 0), 06 <0<0 d - [Ab) 

For 06 < < 4> d , we have sin 2 = 2~ I ^ d ~ 1 ' > = e~ /ln2 /( d_1 \ which shows that there are 
two important cases in terms of the number of vectors. If log M <ti d (A/" 2 d ), a situation 
we refer to as a sparse collection of random vectors, we have I < log M <C d, giving 



9 Jln2 , 7T /Jin 2 , k . 



over the entire range 06 < < <fid- in particular, we have 06 ~ 7r/2 — J\nM/{d — 1). The 
number of groups increases so fast as retreats from ir/2 that for a sparse collection, there is 
a group for each vector when the radius is still quite close to tt/2. In contrast, if log M 3> d 
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(Af ^> 2 d ), which we call a dense collection of vectors, then 0& ~ sin b = 2~ log - /v '/ 2 ( d ~ 1 ) <c 1, 
meaning that to get to one vector per group, the radius 4>b must be small. 

When we turn to the entropy of the groups, it becomes clear that there is yet another 
critical angle, C , the angle at which the number of vectors per group equals the Hilbert- 
space dimension, i.e., A/y(0 c ) = d or J(0 C ) = log Af — log d. For > C , there are sufficiently 
many vectors in each group to explore all the available Hilbert-space dimensions, so the 
entropy is close to the entropy Hd(4>) of a mixture of vectors distributed uniformly within 
a sphere of radius in d dimensions. In contrast, for 0& < < C , the vectors in a group 
can explore roughly only A/y(0) — AA(sin 2 0) d_1 = 2^ I Af dimensions, thus giving an entropy 
close to Hj^ v ^((j)). Finally, for < fe , there is only one vector per group, so H = 0. 

Our main interest is the relation between H and /, so we eliminate the radius from the 
above expressions. The region < 06 gives H = and / = log Af. For 0& < < C , i.e., 
log Af > I > log Af — logd, we have 



H = H MvW (<f)) = H 2 (X) + (1 - A) log(2-W- l) , (A7) 



where 



A = 1 - a-'/*") = 1 - (A8) 



Finally, for C < < 0^, i.e., logA/" — logrf > / > 1, we have H = i?d(0), with 
Summarizing, we have 



Ao = 1 - ( Ll2- 1 '^ . (A9) 



H= fi^ vW (0) = ^ 2 (A) + (l-A)log(2-W-l) , logAf>/>logAf-logd, fA1Q) 
l# d (0) = if2(Ao) + (l-A o )log(d-l), logAf-logd>/>l. 

with A and Ao given by Eqs. (jA8|) and (jA9|) . Equation ()A10|) is the approximate trade-off 
relation we are seeking. 

The important part of the trade-off relation is the part that is independent of the number 
of random vectors, i.e., for 1 < J < log A/" — logrf. Notice that to investigate this region, 
we need Af ^> d, but we do not need Af so large that the random vectors sample generic 
vectors, which would require at least Af ~ 2 d vectors, i.e., a dense collection. We emphasize 
that we do not need a dense collection of vectors to investigate the important part of the 
trade-off relation. 

Before going further, it is useful to put the trade-off relation ()A10|) in other forms, which 
can be easily specialized to the case of a sparse collection of vectors. For the second case, 
which is the case of interest, we can write 

H d (<j>) = \ogd- ^([d(l - 2- 1 '^) + 2- 1 /^} 

x log[d(l - 2- I ' l - d -^) + 2- I l {d - 1 ^} - n- 1 '^) . (All) 

For a sparse collection of vectors, for which I < log Af — \ogd < log AT <C d, or anytime we 
have / <C d, we can approximate this by 

H d (4>) =logd-^((l + /ln2)log(l + /m2) - /) . (A12) 
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We can manipulate the first case in Eq. (jA10|) in a similar way: 

H Mv{4>) {4>) = log AT - I - A log^Qj - (1 - A) log ( j± . (A13) 

The factor 2 / /A/ r increases from 1/d at / = log A/" — logrf to 1 at / = log M . For a sparse 
collection, we can approximate A by 

, 2 1 Jln2 / 2 J \ fk . 

A = A? + d=T( 1 -j5'J " (A14) 

The second term is always small. When the first term dominates, the second two terms in 
Eq. (jA13|) are small. When the first term is as small or smaller than the second, the second 
two terms in Eq. (jA13|) are again small. Thus for a sparse collection, it is always a good 
approximation to use Hj\f v ^)((f)) = log Af — I. 

The conclusion of these considerations is that for sparse collections, the trade-off rela- 
tion (|A10|) is well approximated by 

logA/"-/, logjV > I > logjV- log d, 

logd- -((l + /ln2)log(l + /ln2) - /) , logJV- logd > I > 1. ^ A15 ^ 

This is the form of the trade-off relation that we use in Sec. Ill A "21 When d is large, it is quite 
a good approximation for sparse collections of random vectors, certainly more than adequate 
given the approximate character of the entire sphere-grouping model. These approximate 
expressions are poorest at the knee between the two behaviors, which is also where the 
approximate treatment of the grouping is at its worst. 
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APPENDIX B: ENTROPY OF EQUAL PARTITIONS OF PROJECTIVE 

HILBERT SPACE 

Let \ej), j = 1, . . . , d, be an orthonormal basis for a <i-dimensional Hilbert space, and let 

n 

£f = £l e i>< e ;l (Bl) 
i=i 

be the projector onto the subspace S + spanned by the first n vectors, 

n+m 

P-= E ( B2 ) 

j=n+l 

be the projector onto the subspace £_ spanned by the next m vectors, and P Q = 1 — P + — P_ 
be the projector onto the subspace S spanned by the remaining d — n — m vectors. An 
arbitrary normalized vector can be expanded uniquely as 

|V>) = cosf(cos0|x) + sin%)) + sin£|0) , (B3) 

where \x) £ S+, \t)) G S-, and \<f>) G So are normalized vectors. The angle ^ is the Hilbert- 
space angle between and the span of S + and 5_, and 9 is the Hilbert-space angle between 
the projection of into the span of S+ and i.e., (P + + P-)\ip), and the subspace S + . 



29 



We are interested in the density operator formed from all pure states whose projection 
into the span of S + and S- is closer to S + than an angle 6, 

p = JV\ d5 M -i|^)^| , (B4) 

where JV is a normalization factor. Here and throughout dSj denotes the standard inte- 
gration measure on the j-sphere, and Sj = J dSj is the volume of the j-sphere. This region 
of states is the analogue of the intersection in three real dimensions of a wedge of opening 
angle 2G with the unit sphere. It is clear that p is invariant under unitary transformations 
that are block-diagonal in the three subspaces, which implies that p has the form 

p = A+P+ + A_P_ + A P • (B5) 

Our job is to determine the three eigenvalues, A± and A , which satisfy 

nX + + mA_ + (d — n — m)\ = 1 . (B6) 

It turns out that A = l/d, as we show below, so we have 

X- = -Jl + -(l-d\ + )) . (B7) 



d V m 

A small change in \ip) can be written as 

|d^) = d£(— sin£(cos0|x) + sin^|?y)) + cos£|0)) + sin£|d0) 

+ cos£(d0(-sin0|x) + cos%)) + cos6>|dx) + sin^dr?)) . (B8) 

This gives a line element on normalized vectors, 

ds 2 = (<ty|<ty) = d£ 2 + sin 2 £(d0|d0) + cos 2 e(d^ 2 + cos 2 #(dx|dx> + sin^d^d??)) , (B9) 

and a corresponding volume element on the (2d — l)-sphere of normalized vectors, 

d<S 2d _i = sin 2(d - n - m) - 1 ^cos 2(n+m) - 1 ^d^ 

x cos 2 "- 1 ^ sin 2 ™"^ d6 d«S 2{d _ n _ m) _ 1 d5 2n _! dS 2m - 1 . (BIO) 

Normalizing the density operator gives 

1 = tr(p) = JV ( dS M _i 
Je<e 

x [°d6 cos 2n - 1 0sin 2m - 1 . (Bll) 
Jo 

We first verify that X = l/d. Letting |e ) be any normalized vector in So, we have 
A = (e |p|e ) = JV / d5 2d _i|(e |V')| 2 

Jo 

x J\e cos 2 "-^sin 2m -^ J dS 2{d - n -m)-i \(e \<P)\ 2 • (B12) 
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Using 

/ dS 2{d ^ m) ^ |(e o |0)| 2 = ^ (B13) 

J a — n — m 

and the expression for the normalization constant from Eq. (|B11|) and changing integration 
variable to u = sin 2 £, we get 



i [ duu d - n - m (l-u 

\ _ 1 £0 

~~ YZ. — Z 7 1 

CZ Tl Tfl I i d— n — m —i/-. \ n _|_ m _i 



n+m— 1 



dllU m yl — u) 

1 T(d-n-m + l)T(n + m)/T(d + 1) 
d — n — m T(d — n — m)Y{n + m)/T(d) 
1 

d ' 



(B14) 



Similarly, to find A + , we let \e + ) be any normalized vector in S + and write 
A + = (e + |p|e + ) =jr\ d^Ke+l^)! 2 

= ^«S 2(d _ n _ m) _ 1 «S 2m _ 1 ^ 7r/2 d^ sin^-^-^cos 2 ^^ 

x J\e cos 2n+1 0sui 2m - l J dS 2n -i |(e+|x)| 2 • (B15) 



Using 



dS 2n - X \{e+\ X )\ 2 = — (B16) 
n 

and the expression for the normalization constant and changing integration variables to 
u = sin 2 £ and v = sin 2 ^, we get 



A. 



1 rsin 2 e 

i / d«« d_n-m - 1 (l -u) n+m / dvv m -\l-v) n 
1 Jo Jo 

" ^ duu d - n - m - l {l-u) n+m - 1 [ Sm2e dw^ll-v) 11 ' 1 







lT(d-n-m)T(n + m + l)/T(d+l) J Q 



sin 2 

dv v m ~\i - v y 



n r(d-n-m)T(n + m)/T(d) ^ ° dvv^ 1 {I - v) n ~ l 



n + m 



o 



sin 2 

6.vv m -\\-v) n 



nrl /-sin 2 © 

dvv m -\l-v) n - 1 



(B17) 



We now specialize to the case of interest, n = m and = 7r/4, so that p is constructed 
from pure states occupying one of two halves of Hilbert space: 

/■1/2 

« / du^-^l-u)" 

V = J% • (B18) 

; dw^Hl-vY' 1 
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The integrals can be evaluated as 

Jo V ; 2(2n)! V / 

/ d^^a-^r- 1 = m — d_ . (Big) 

jo 2(2n — 1)! 

Plugging these results into Eqs. (jB18|) and (jB7|) . we get 

1 / r(n + l/2) \ _ 1 / (2n)! \ 

When n = 1, we get A + = 3/ 2d and A_ = l/2d, and when n — 2, X + — ll/8d and A_ = 5/8d. 
For large n (and d), we can use Stirling's formula to write 

A± « i (l ± -Ll . (B21) 



The von Neumann entropy of p can be put in the form 

H = — n\ + log A + — nA_ log A_ — (d — 2n)\ log A: 



o 



9 n 

= \ogd-—{l-H 2 (dX + /2)) . (B22) 

For fixed d, this is a decreasing function of n. For large n (and cf), we can use Eq. (IB21)) to 
write H 2 {d\ + /2) w 1 - l/2?rn In 2 and 

1 46 

if « logrf - — — = logd - — - . (B23) 
iraln 2 d 
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